( 


AD- A 143  528 


XFOSR-TR.  8  4  ..  o  5  g  o 


CAR-TR-5  3 
CS-TR-1 33  G 


F49620-83-C-0032 
March  1934 


DEGREE  OF  ADJACENCY 

or  surrounded:: ess 

Azricl  Rosenfeld 

Center  for  Automation  Research 
University  of  Maryland 
College  Park,  MD  20742 

Reinhard  Klette 

Image  Processing  Laboratory 
Friedrich- Schiller  Univors.it  a  c 
DDR-69  Jena 

German  Democratic  Republic 


COMPUTER  VISION  LABORATORY 


CENTER  FOR  AUTOMATION  RESEARCH 


UNIVERSITY  OF  MARYLAND 

COLLEGE  PARK,  MARYLAND 
20742 


8  4  07  24  049 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fHTien  Data  Entered) 

REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

■itfllfliii 

3.  RECIPIENT'S  CATALOG  NUMBER 

S _ 

4.  TITLE  (and  Subtitle) 

DEGREES  OF  ADJACENCY  OR  SURROUNDEDNESS 

5.  TYPE  OF  REPORT  4  PERlOO  COVERED 

Technical 

6.  PERFORMING  ORG.  REPORT  NUMBER 

CAR-TR-5  3:  CS-TR-13a_0  _ 

7.  AUTHORf*) 

Azriel  Rosenfeld 

Reinhard  Klette 

8.  CONTRACT  OR  GRANT  NUMBERf«J 

F49620-83-C-0082 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Center  for  Automation  Research 

University  of  Maryland 

College  Park,  MD  20742 

to.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  &  WORK  UNIT  NUMBERS 

OtlCSjfr 

It.  CONTROLLING  OFFICE  NAME  AND  AODRESS 

Math.  &  Info.  Sciences,  AFOSR/NM 

Bolling  AFB 

Washington,  DC  20332 

12.  report  DATE  / 

March  1984/ 

13.  NUMBER  OF  PAGES 

32 

14.  MONITORING  AGENCY  NAME  &  ADDRESS^//  different  from  Controlling  Office) 

15.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

IS  a.  OECLASSIFI  CATION  DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


Approved  for  public  release?  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  abatract  entered  in  Block  20,  if  different  from  Report) 


18.  supplementary  notes 


19.  KEY  WORDS  (Continue  on  reverse  aide  if  necessary  and  identify  by  block  number) 

Geometry  Adjacency 

Digital  geometry  Surroundedness 

Fuzzy  geometry 
Spatial  relations 


20.  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

Definitions  of  the  degree  of  adjacency  of  two  regions  in  the 
plane,  and  the  degree  of  surroundedness  of  one  region  by  another, 
are  proposed.  Some  elementary  properties  of  these  concepts  are 
established,  and  it  is  also  shown  that  they  have  natural  generali¬ 
zations  to  fuzzy  subsets  of  the  plane.  Applications  of  the  proposed 
measures  to  digital  polygons  are  demonstrated  and  fast  algorithms 
for  computing  these  measures  are  given. 


DD 


FORM 
1  JAN  73 


1473 


EDITION  OF  1  NOV  «S  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  This  page  fl»>.en  Data  Entttad) 


CAR-TR-53 

CS-TR-1380 


F49620-83-C-0082 
March  1984 

DEGREE  OF  ADJACENCY 
OR  SURROUNDEDNESS 

Azriel  Rosenfeld 

Center  for  Automation  Research 
University  of  Maryland 
College  Park,  MD  20742 

Reinhard  Klette 

Image  Processing  Laboratory 
Friedrich-Schiller  Universitat 

DDR-69  Jena  ,  • 

German  Democratic  Republic  ^  ^ 


ABSTRACT 
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by  another,  are  proposed.  Some  elementary  properties  of 
these  concepts  are  established,  and  it  is  also  shown  that  they 
have  natural  generalizations  to  fuzzy  subsets  of  the  plane. 
Applications  of  the  proposed  measures  to  digital  polygons  are 
demonstrated  and  fast  algorithms  for  computing  these  measures 
are  given. 
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1.  Introduction 


In  describing  a  picture,  one  often  needs  to  specify  geo¬ 
metric  relations  among  the  regions  of  which  the  picture  is 
composed.  For  a  review  of  such  relations  and  their  measure¬ 
ment  in  digital  pictures,  see  [1]. 

Adjacency  is  an  important  relation  between  regions.  In 
a  digital  picture,  sets  S  and  T  are  adjacent  if  seme  border  pixel  of 
S  is  a  neighbor  of  some  border  pixel  of  T;  in  the  Euclidean  plane, 
regions  S  and  T  are  adjacent  if  their  borders  intersect. 

Note  that  this  relation  is  not  quantitative;  S  and  T  are  not 
considered  adjacent  even  if  they  are  very  close  to  one  another, 
and  it  also  makes  no  difference  whether  they  are  adjacent  at 
one  point  or  at  many  points.  In  Section  2  we  propose  a  quan¬ 
titative  definition  of  adjacency  which  does  take  these  factors 
into  account. 

Another  important  region  relation  is  surroundedne ss .  We 
assume  that  all  pictures  are  of  finite  size;  the  region  of 
the  plane  outside  a  picture  is  called  the  background.  We 
say  that  S  surrounds  T  if  any  path  from  T  to  the  background 
must  intersect  S.  This  definition  too  is  nonquantitative . 

In  Section  3  we  propose  two  ways  of  defining  the  degree  to 
which  S  surrounds  T. 


2 .  Quantitative  adjacenc 


In  Section  2.1  we  define  quantitative  adjacency  for 
regions  in  the  Euclidean  plane.  In  Section  2.2  we  general¬ 
ize  this  definition  to  fuzzy  subsets  of  the  plane,  and  in 
Section  2.3  we  discuss  quantitative  adjacency  of  subsets  of 
a  digital  picture. 

2.1.  Euclidean  regions 

Let  CQ  be  a  rectifiable  simple  closed  curve  in  the  plane, 

and  let  C^,...,C  be  rectifiable  simple  closed  curves  not 

crossing  one  another  and  contained  in  the  interior  of  CQ . 

According  to  the  orientation  of  CQ,  the  closed  set  CgUCg  will 

be  either  a  bounded  set,  or  the  infinite  plane  except  for  a 

bounded  set  (infinite  case).  For  C^,...,Cn  we  assume  that  Clf 

—  ,C_  are  bounded  sets.  Then  (C.UCn)  -  (C-.U  .  .  .Uc  )  is  called 
n  0  U  1  n 

a  region;  C„  is  called  its  outer  border,  and  C, , . . . ,C  are 
— -  u  -  1  n 

called  its  hole  borders.  Note  that  in  the  infinite  case, 

there  is  no  unique  distinction  between  the  outer  border  and 

the  hole  borders  because  the  outer  border  may  be  considered 

to  be  a  hole  border  itself.  The  perimeter  of  the  region  is  the 
n 

sum  Z  | C . |  of  the  lengths  of  its  borders. 
i=0  1 

Intuitively,  two  regions  S  and  T  are  (somewhat)  adjacent 
if  some  border  of  S  (nearly)  touches  some  border  of  T;  the 
degree  of  adjacency  depends  on  how  nearly  they  touch  and 
along  how  much  of  their  lengths  they  do  so.  The  borders 
nearly  touch  if  they  are  close  to  one  another,  as  illustrated 
in  Figure  la-b;  note  that  S  and  T  are  allowed  to  overlap. 


/; 


Note,  however,  that  not  all  cases  in  which  borders  are  close 
to  each  other  imply  near-adjacency,  as  shown  in  Figure  lc. 

The  difference  is  that  in  Figures  la-b,  the  shortest  paths 
between  the  close  borders  lie  outside  both  regions  or  inside 
both  of  them,  while  in  Figure  lc  these  paths  lie  inside  one 
region  and  outside  the  other.  It  also  seems  plausible  that 
only  line-of  sight  paths  should  be  counted  in  defining 
adjacency;  in  Figure  2a,  the  left-hand  .edge  of  T  should  con¬ 
tribute  to  its  degree  of  adjacency  to  S,  but  its  other  edges 
should  not,  and  similarly  in  Figure  2b,  the  parts  of  the  bor¬ 
der  of  the  concavity  in  T  from  which  S  is  not  visible  should 
not  contribute.  Finally,  note  that  quantitative  adjacency 
is  not  symmetric;  in  Figure  3(a-b),  S  is  highly  adjacent  to 
T,  since  much  (or  all)  of  its  border  nearly  coincides  with 
the  border  of  T,  but  T  is  not  as  highly  adjacent  to  S,  since 
only  a  small  fraction  of  its  border  coincides  with  that  of  S. 

Based  on  these  considerations,  we  define  the  degree  of 
adjacency  of  S  to  T  as  follows:  Let  P,Q  be  any  border  points 
of  S  and  T,  respectively.  If  P^Q,  we  say  that  the  line  seg¬ 
ment  PQ  is  admissible  ("with  respect  to  (S,T)"  understood)  if 
its  interior  lies  entirely  outside  both  S  and  T  or  entirely 
inside  both  of  them.  If  P=Q,  we  call  PQ  admissible  if  the 
(signed)  normals  to  the  borders  of  S  and  T  at  P  do  not  point 
in  the  same  direction.  Let  dp  be  the  length  of  the  shortest 
admissible  line  segment  PQ  having  P  as  an  endpoint;  if  no 
such  segment  exists,  let  d  =°°.  Then  we  define  a  (S,T)=  -y 
where  the  integration  is  over  the  border  9S  of  S.  For  example 


if  S  and  T  are  two  squares  of  size  a*a  with  distance  b  between 
them,  then  a  (S  ,T)  =a  (T ,S  )  =  ^  .  Or,  if  S  is  a  square  of  size 
a*a  located  at  the  center  of  a  square  T  of  size  bxb,  defining 
an  infinite  region  T,  then  a(S,T)=  (j^a)/2  '  SO  that  a^s'T^  "*■ 

4b  for  a  -*•  b. 

According  to  this  definition,  a  border  point  P  of  S  con¬ 
tributes  maximally  to  a(S,T)  if  it  also  lies  on  the  border 
of  T  (and  the  conditions  for  the  case  P=Q  are  met) ,  since  in 

this  case  d  =0  and  ,  ,  =1;  and  it  does  not  contribute  at  all 

P  dp+l 

if  no  admissible  segment  PQ  exists  (e.g.,  if  the  border  of  T 
is  not  visible  from  P)  ,  since  in  this  case  d  =°°  and  3-1--  =0 

p  dp+l 

Since  dp20,  in  any  case  we  have  t^jSI;  thus  a(S,T)£/l  dP= 

aP  8S 

p(S),  the  perimeter  of  S.  If  desired,  one  can  normalize 
a  (S ,T)  by  dividing  it  by  p(S);  it  then  lies  between  0  (not 
at  all  adjacent)  and  1  (maximally  adjacent).  For  example, 
a  hole  in  a  region  is  maximally  adjacent  to  that  region;  see 
Proposition  2.2.  We  could  have  used  a  different  function 
f  (dp) ,  in  place  of  ^  ,  in  defining  a(S,T);  the  essential 

requirements  are  that  f  be  a  monotonically  decreasing  func¬ 
tion  of  dp,  and  that  f(0)=l,  f(“)=0. 

Proposition  2.1.  For  regions  S  and  T,  a(S,T)=0  if  and  only 
if  Sct. 

Proof ;  If  S£T  holds,  there  are  no  admissible 
segments  PQ.  Note  that  where  the  borders  touch,  the  signed 
normals  of  S  and  T  point  in  the  same  direction.  Conversely, 


if  S£T,  then  by  the  definition  of  a  region,  there  must 
exist  a  border  arc  of  S  at  every  point  of  which  there  is 
an  admissible  segment,  so  that  a(S,T)^0.  || 

Proposition  2.2.  For  regions  S  and  T,  a(S,T)=p(S)  if  and 
only  if  either  S  is  bounded,  lies  inside  a  hole  in  T,  and 
its  border  is  identical  to  the  border  of  that  hole;  or  S 
is  unbounded,  T  lies  inside  a  hole  in  S,  and  the  border  of 
S  is  identical  to  the  outer  border  of  T. 

Proof ;  If  these  conditions  hold,  we  have  dp=0  at  every 
border  point  of  S.  Conversely,  under  no  other  circumstances 
can  the  entire  border  of  S  coincide  with  a  part  of  the  border 
of  T.  Note  that  a  region  cannot  consist  of  several  isolated 
parts;  thus  it  must  be  contained  in  one  hole  only.  || 

The  definition  of  a(S,T)  can  easily  be  extended  to  vari¬ 
ous  types  of  sets  other  than  regions.  For  example,  for  a 

-  1 

single  point  P  we  can  define  a(P,T)  as  0  if  P€T,  and  ^  fp  t) +1 
if  P/f,  where  d(P,T)  is  the  distance  from  P  to  T.*  (It  follows 
from  this  definition  that  a(P,T)=l  when  P  is  on  the  border  of 
T.)  Conversely,  for  a  single  point  Q  we  can  define  a(S,Q) 
using  our  original  definition  for  sets  S  and  {Q/;  note  that 
here  too  there  are  no  admissible  segments  if  QiS ,  but  that 
otherwise  a(S,Q)  is  obtained  by  integration  over  the  part 
of  the  border  of  S  visible  from  Q.  Similarly,  we  can  define 

*-his  definition  is  not  exactly  analogous  to  the  one  for  regions 
a  single  point  has  zero  border  length,  so  that  integrating  over 
it  should  always  give  zero.  The  analogy  would  be  better  if,  in 
the  region  definition,  we  normalized  a(S,T)  by  dividing  p(S) 


a(S,T)  if  S  or  T  is  an  arc;  the  details  are  left  to  the  reader. 

If  a(S.,T .)  are  all  defined  for  S  ,...,S  and  T.  ,  .  .  .  ,T  ,  we 
3  m  n  1  m  1  n' 

can  also  define  a(  U  S.  ,  U  T  .  )  ;  again,  the  details  are  emitted. 

i=l  j=l  3 

Proposition  2.3.  If  T'ct  and  SflT=0,  then  a  (S ,  T  '  )  ia  ( S  ,T )  . 


Proof :  Let  PQ  be  any  admissible  segment  in  the  definition 

of  a(S,T').  Let  R  be  the  first  point  in  which  PQ  meets  T; 
then  PR  or  a  shorter  line  segment  is  admissible  in  the  defi¬ 
nition  of  a(S,T) .  (Note  that  R^P  since  sHt=0 . )  If  dp,dp  are 
the  lengths  of  the  shortest  admissible  segments  in  the  defi¬ 
nitions  of  a(S,T)  and  a(S,T'),  respectively,  we  thus  have 
dp-d^.  (Note  that  for  some  P's  there  may  be  admissible  seg¬ 
ments  with  respect  to  T  but  not  with  respect  to  T ' . )  Thus 


/; 


as  dp+1 


dP  s  j 


as  dp+1 


dP. 


Proposition  2.4.  If  P,P'/T  and  d  (P ' ,T) ^d  (P ,T)  ,  then  a(P',T)  2 
a(P,T),  but  it  is  not  necessarily  true  that  a (T , P ' ) 2a  (T , P)  . 


Proof:  Evidently  a(P',T)  -  d(p,  ;t)+T  2  dTF^T+1  =  a(P'T)* 

On  the  other  hand,  let  P,  P',  and  T  be  as  shown  in  Figure  4; 
then  evidently  (for  P  sufficiently  close  to  P')  the  contribu¬ 
tions  of  side  b  to  a(T,P)  and  a(T,P')  are  approximately  equal, 
but  sides  a  and  c  do  not  contribute  to  a(T,P'),  so  that  its 
total  contribution  is  smaller  (see  example  2.2). 


Proposition  2.5.  For  P/T,  a(P,T)  is  a  continuous  function  of 
the  position  of  P,  but  a(T,P)  is  not. 

Proof:  See  the  proof  of  Proposition  2.4  and  the  example 


given  there. 


2.2.  Fuzzy  subsets 

Our  definition  of  a(S,T)  can  be  generalized  to  the  case 
where  S  and  T  (call  them  y  and  v)  are  bounded  fuzzy  subsets  - 
that  is,  functions  defined  on  the  plane,  with  values  in  [0,1], 
and  equal  to  zero  outside  a  bounded  region  B.  The  desirability 
of  defining  geometric  concepts  in  the  fuzzy  case,  so  that  they 
can  be  measured  without  having  to  first  crisply  segment  a  pic¬ 
ture,  is  discussed  in  [ 2 ] _ 

We  assume  in  what  follows  that  u  (and  v*  .xe  "piecewise 
constant"  in  the  following  sense:  We  part  ^on  B  into  a 
finite  number  of  regions  whose  interiors  a  ^  is joint,  and 
where  the  border  of  each  region  is  contained  in  the  union  of 
the  borders  of  the  other  regions.  Let  6  be  the  union  of  all 
borders  of  these  regions.  In  the  interior  of  each  region, 
y  has  constant  value;  at  each  point  of  a  border,  it  has 
one  of  the  neighboring  interior  values.  [Another  case 
of  interest  is  that  in  which  y  and  v  are  "smooth,"  i.e., 
everywhere  differentiable;  note  that  we  can  approximate  a 
piecewise  constant  u  by  a  smooth  v  which  is  constant  except 
near  the  borders,  where  it  changes  rapidly  from  one  constant 
value  to  another.  "Smooth"  versions  of  the  definitions  in 
this  section  could  be  given,  using  derivatives  in  place  of 
differences . ] 

If  P/Q,  we  call  the  segment  PQ  admissible  ("with  respect 


to  (.  - ,  v )  "  understood)  if: 


a)  P  is  on  a  border  of  B,  and  Q  on  a  bo^cr  cf  g. 

We  assume  that  or 1 y  L«o  of  the  constant  regions 
r-z  u(v)  meet  at  P(Q).  (More  than  two  regions 
met  only  at  a  finite  number  of  points,  and  these 

can  be  ignored  in  defining  degree  of  adjacency.)  ( 

b)  Let  R^P  be  a  point  of  PQ  such  that  u  changes 
value  at  R  as  we  move  from  P  to  Q  (there  can  only 

be  finitely  many  such  R's) ,  and  let  AR  be  this  ^ 

change  in  value.  Let  the  values  of  g  at  the  two 
regions  that  meet  at  P  be  a  and  b,  where  b  is  the 

value  on  PQ  near  P,  and  let  Ap=a-b.  We  assume  ^ 

that  Ap  and  all  the  AR's  have  the  same  sign,  and 
that  | Ap | > | Ar |  for  all  R. 

c)  Let  vR  and  be  defined  analogously  to  iR  and  lp  4 

in  (b)  ,  with  j  replacing  t  and  the  roles  of  P  and 

Q  reversed.  We  assume  that  7  and  all  the  7R's  have 

the  same  sign;  that  this  is  the  same  sign  as  in  (b) ;  4 

and  that  '  7  for  all  R. 

y  k 

Conditions  (b-c)  mean  that  the  changes  in  ..  as  we  move  from 

P  to  Q  are  all  in  the  same  direction,  and  the  change  "at  P"  is  4 

the  largest  of  them;  and  analogously  for  the  changes  in  as 
we  move  from  Q  to  P .  Intuitively,  this  means  that  the  border 

points  P,Q  of  6  are  facing  toward  or  away  from  each  other  4 

(since  the  changes  have  the  same  signs  in  both  cases)  ,  and 

that  no  "stronger”  border  points  (at  which  larger  changes 

occur'  lie  between  them,  so  that  they  are  within  "line  of  sight.”  * 


i 


It  is  easily  seen  that  if  y  and  v  are  crisp,  these  conditions 
reduce  to  the  definition  of  admissibility  given  in  Section 
2.1  for  P^Q. 

When  P=Q,  we  call  PQ  admissible  if  P  is  a  border  point 
of  3,  and  the  changes  in  y  and  v  at  P  (in  a  fixed  direction 
from  one  region  of  the  partition  of  B  to  the  other  touching 
it  at  P)  have  opposite  sign.  In  this  case,  let  Ap  and  Vq 
be  the  changes  in  the  values  of  y  and  v  at  P,  defined  as 
in  the  preceding  paragraph,  where  Ap'7^-0  is  assumed. 

For  each  P,  let  g(P)  =  _ SU^  AP  7Q  ,  where  d(P,Q) 

PQ  admissible  d(P,Q)+l 

is  the  distance  from  P  to  Q.  Note  that  the  numerator  is  always 
postiive,  since  the  changes  in  u  and  v  at  P  and  Q  have  the  same 
sign;  and  that  |Ap|  and  |Vg]  each  ^1,  so  that  the  numerator  is 
in  the  interval  (0,1].  Evidently,  in  the  crisp  case  the  numera¬ 
tor  must  be  1,  and  the  sup  is  achieved  when  the  denominator  is 
as  small  as  possible,  so  that  g  (P)  is  the  same  as  ^ j-  of 
Section  2.1.  It  is  understood  that  g(P)=0  if  no  admissible  PQ 
exists . 

Note  that  our  definition  of  g(P)  involves  a  tradeoff  between 
the  border  strengths  (=sizes  of  changes)  at  P  and  Q  and  the  dis¬ 
tance  d(P,Q);  the  sup  may  arise  from  weak  changes  that  are  close 
together  (or  even  coincide) ,  or  from  stronger  changes  that  are 
farther  apart.  The  nature  of  the  tradeoff  can  be  manipulated 
by  using  some  other  monotonic  function  of  d(P,Q)  in  place  of 
^  ;  compare  the  remark  about  f  (P)  in  Section  2.1. 


Finally,  we  define  a(y,v)=  / g(P)dP,  where  the  integra- 

B 

tion  is  along  all  borders  6  of  the  partition  of  B.  We 
leave  it  to  the  mcerested  reader  to  consider  extensions 
of  this  definition  (e.g.  ,  to  define  a(P,v)  and  a(y,Q), 
where  P  and  Q  are  points),  and  to  investigate  the  possi¬ 
bility  of  fuzzy  generalizations  of  the  Propositions  in 
Section  2.1.  We  prove  here  only 


Proposition  2.6.  a  (y  ,  v)  sp  (y )  ,  the  (fuzzy)  perimeter  of  y. 
Proof:  We  recall  [3]  that  P(u)  is  just  the  sum  of  the  lengths 
of  the  border  arcs  of  3  at  which  pairs  of  regions  of  y  meet, 
each  multiplied  by  the  absolute  difference  in  value  between 
that  pair  of  regions.  Now  this  difference,  at  a  given  border 
point  P,  is  just  Ap.  Thus 


P  (y)  =  /  |  A  |  dP  >  /  g  (P )  dP 

3  P  B 


for  any  v,  since  in  any  case  g(P) 
point  Q  with  |VQ|<1  and 


=  dTPToT+T  for  a  certain 
s 1,  while  Ap-VQ=|Apl  *|vq 


since  they  have  the  same  sign.  || 


The  fuzzy  generalization  of  Proposition  2.3  is  false; 


yAv-0  and  v'sv  does  not  imply  a  (y  ,  v ' ) ia  (y , v ) .  In  Figure  5, 

PQ  is  admissible  for  V,  and  g  (P )  =  -p—— j-  .  With  respect  to  v, 

however,  the  steps  in  value  are  all  1/n,  except  for  the  last 
2 

which  is  — ;  the  only  possible  maximal  values  of  g(P)  are  thus 

d+T  anc*  d+n-1  ‘  ^  n>2  anc*  d+n-l<  (d+1 )  n ,  these  values  are 

both  smaller  than  that  for  v'.  This  counterexample  would  break 


down  if  we  used  a  different  definition  of  g(P),  such  as 

Ap  t  Ap .  ^  y 

^ d  ( P , Q) 1 1  +  R?PQ  ^d (R,Q) +1 ^  ^  ^d (Q , P ) +1  +  R£Qp  ^d (R , P ) +1 ^  ^  ; 


but 


we  will  not  pursue  this  alternative  approach  here. 


Subsets  of  digital  pictures  may  be  considered  from 

different  points  of  view,  e.g.,  as  sets  of  grid  points, 

as  sets  of  cells,  or  as  digital  polygons.  For  purposes  of 

defining  quantitative  adjacency  in  the  digital  (crisp) 

case,  it  is  convenient  to  deal  with  digital  polygons. 

In  a  digital  simple  polygon  S  =  <Po ' P1 ' ' ' * ' Pn> '  for 

k  =  0,1,..., n,  the  PR  are  all  grid  points  with  integer 

coordinates;  points  and  pk+1  are  8 --neighbors  (for 

p  =p  ) ;  in  relation  to  the  interior  of  S  the  sequence 
n+1  0 

Po'Pl'...'Pn  haS  clockwise  orientation;  and  the  border  of 
S  is  non-crossing,  i.e.,  S  is  a  simple  polygon  in  the 
usual  sense.  Because  of  the  fixed  orientation,  finite 
and  infinite  digital  simple  polygons  can  both  be  defined 
in  this  way.  For  S  =  <Pq  ,  ,  •  •  •  ,  pn>  >  complementary 

polygon  S  is  given  by  <Pn,pn_1 , • • •  ,p0> * 

Because  digital  simple  polygons  are  regions  as  defined 
in  Section  2.1,  the  degree  of  adjacency  a(S,T)  is  defined 
for  digital  polygons  S  and  T.  However,  for  the  needs  of 
picture  processing  or  computer  graphics,  a  more  specifically 
digital  approach  will  be  used.  For  this  purpose,  the  sat 
of  border  points  BP(S)  of  a  digital  polygon  S  will  be  re¬ 
stricted  to  grid  points  on  the  (real)  border  of  S,  i.e., 

BP  (S)  =  v  PQ  ,  pL ,  •  •  •  '  pn  ^  for  s  =  <P0'P1"  •  ‘  ,Pn>'  N°W'  admis~ 
sible  line  segments  are  defined  as  in  Section  2.1  for 
PtBP(S)  and  QtBP (T) ,  for  the  digital  polygons  S  and  T.  (At 


a  vertex  of  a  digital  polygon,  the  (signed)  normal  is  defined 
to  be  the  bisector  of  the  vertex  angle.) .  Let  AL(S,T)  be  the 
set  of  all  admissible  line  segments,  from  border  points  of  S 
to  border  points  of  T.  Note  that  AL (S ,T) ^AL (T , S)  for  almost 
all  digital  polygons  S  and  T;  AL (S , S ) =BP  (S) ,  and  AL(S,S)=0 
(the  empty  set) .  For  digital  polygons  S  and  T,  the  digital 
degree  of  adjacency  is  defined  to  be 

/  2  1/  (l+dp)  ,  if  AL(S,T)^0 

I  P€BP (S) 


WS'T>  ■ 


0  otherwise. 


Note  that  all  the  properties  given  in  Section  2.1  for  function 
a  are  true  for  a^  too.  The  property  a(S,T)£p(S),  the  peri¬ 
meter  of  S,  is  replaced  by  adig(S,T)  -  card  BP(S),  which  may  be 
considered  to  be  the  digital  perimeter  of  S. 

Proposition  2.7.  For  digital  polygons  S,T,  a^  (S  ,T)  =card  BP(S) 
if  and  only  if  T=S. 


Proof :  Because  of  AL (S , S) =BP (S)  we  have  dp=0  for  all  PtBP(S), 

Thus,  a,.  (S,S)=card  BP(S).  On  the  other  hand,  if  T^S  then 
dag 

there  exists  at  least  one  point  P^BP(S)  with  dp>0 . 

The  normalized  degree  of  adjacency  is  defined  by 

adig(S,T)  =  adig  (s'T>/card  BP(S), 
for  digital  polygons  S  and  T.  The  behavior  of  this  concept  of 


degree  of  adjacency  is  illustrated  by  the  following  two  examples. 


Example  2.1.  Assume  two  convex  digital  polygons  S  and 
with  distance  n  between  them,  as  shown  in  Figure  6.  For 
different  values  of  n,  polygon  Tn  changes  its  position  in 
relation  to  S.  For  example,  for  n=0 ,  S  is  in  centralized  posi¬ 
tion  within  Tq,  and  for  n=7  and  n= -7 ,  S  and  Tn  are  in  touch¬ 
ing  positions.  By  symmetry,  for  n=0,l,2,...  we  have 


adig(S'Tn>'adig(S'T-nl  and  adig<Tn'S)"adig(T-n'S)-  For 
the  normalized  degrees  of  adjacency  we  use  card  BP(S)=16  and 


card  BP  (T  )=24.  For  n=0 , 1 , 2 , . . .  ,  10  ,  the  values  of  a^ .  (S,T  ) 
n  dig  n 

and  ajjig(Tn'S)  are  given  in  Table  1.  These  values  are  gra¬ 
phically  illustrated  in  Figure  7.  As  seen  in  this  Figure, 


there  is  a  somewhat  unbalanced  behavior  of  the  proposed  mea¬ 


sure  a*^  for  intersecting  positions  (-65n^6)  and  non-inter¬ 
secting  positions  (|n|28)  of  the  two  polygons.  Even  for 
the  "most  adjacent"  positions  (n=7)  we  don’t  have  the  maximum 
value,  which  is  realized  for  ] n  j  =2 .  This  behavior  is  due  to 
the  influence  of  border  points  for  which  P=Q,  P€BP(S)  and 
Q €BP  (T^)  .  We  may  change  our  definition  of  a^  by  requiring 
dp=0  if  and  only  if  there  are  border  segments  PP1  and  QQ’ 
of  S  and  T  ,  respectively,  such  that  P^P’  (non-trivial  seg¬ 
ment),  PP '=QQ',  and  the  signed  normals  of  S  and  T^  on  this 
common  border  segment  point  in  exactly  inverse  directions.  In 


all  other  cases  of  P=Q,  P€BP(S) 

The  resulting  modified  function 

_  + 
suits  for  a  ,  • 
dig 


and  Q  £  BP  (T  )  ,  we  let  dp=<» 


a*^^  is  denoted  by 


a  ,  .  . 
dig 


Re- 


are  also  shown  in  Table  1  for  -lOtnrlO. 


Remark.  Convex  digital  polygons  are  restricted  to  be  oc¬ 
tagons  at  most,  where  "convexity"  is  understood  as  in  the 
real  plane.  Then,  for  two  convex  digital  polygons  S  and 
T  with  N=card  BP(S)  +  card  BP(T),  there  exists  an  0(N) 
worst  case  time  algorithm  for  computing  a^g(S,T)  or 
a*^g(S,T).  The  basic  ideas  of  the  algorithm  are  as  follows: 

(i)  Determine  the  upper  and  lower  tangents  on  sUt  dividing 
the  polygonal  border  of  S  (or  T)  into  a  connected  part, 
where  no  admissible  line  segments  to  T  (or  to  S)  are  possi¬ 
ble,  and  a  second  connected  part  containing  all  points  which 
may  contribute  to  ad^g(S,T)  (or  to  a^  (T,S)). 

(ii)  Perform  two  search  procedures,  one  top-down  and  one 
bottom-up,  to  compute  candidate  values  for  dp ,  for  all 
points  P  in  the  interesting  part  of  the  border  of  S,  where 
only  points  in  the  interesting  part  of  the  border  of  T 

have  to  be  considered.  In  both  search  procedures  the  connec¬ 
tion  line  between  two  points  P6BP(S)  and  Q6BP(T)  which  are 
under  consideration  for  computing  dp  moves  monotonically  down 
(or  up)  in  the  interesting  parts  of  the  borders  of  S  and  T. 
During  these  search  procedures  at  most  two  crossings  of 
the  borders  of  S  and  T  are  possible. 

(iii)  For  all  points  P  in  the  interesting  part  of  the  border 
of  S  take  the  minimum  of  both  candidate  values  found  in  (ii) 
to  compute  dp. 

(iv)  Determine  a^g(S,T)  by  using  the  values  of  dp. 

Using  small  examples,  it  can  already  ben  seen  that  in 

fact  two  search  procedures  are  necessary  in  step  (ii)  which 
cooperate  to  give  the  final  result  in  step  (iv) . 


3 .  Quantitative  surroundedness 

In  Sections  3.1  and  3.2  we  present  two  quantitative 
definitions  of  surroundedness  in  the  Euclidean  plane, 
and  also  indicate  how  each  of  them  generalizes  to  the  fuzzy 
case.  Section  3.3  discusses  quantitative  surroundedness 
in  digital  pictures. 

3.1.  Visual  surroundedness 

Let  P  be  a  point  and  T  a  (bounded)  set.  Let  r„(P,T)=l  if 
the  ray  emanating  from  P  in  direction  9  meets  T,  and  rg(P,T)  =  0 
otherwise.  We  define  the  degree  of  visual  surroundedness  of  P 
by  T  as  v(P,T)=~r  j  ^  Tq  (P,T)  dQ  .  (This  integral  might  not  be 
defined  if  T  is  sonetypeof  pathological  set;  but  it  is  evidently 
defined  for  various  types  of  well-behaved  sets,  such  as  regions 
and  arcs. ) 

If  S  is  a  (well-behaved)  set,  we  define  v(s,T)  as 

min  v(P,T).  (Another  possibility  would  be  to  take  the 
P<ES 

"average"  value  of  v(P,T)  for  all  P?S. )  It  follows  that 

v  (S ,T)  is  defined  by  a  border  point  of  S,  i.e.,  v(S,T)  = 

min  v(P,T).  For  the  sets  given  in  Figure  9,  we  have 
P^dS 

v (S,T)=v (P,T)=tan  ^  (5/6 ) /n  =  0 . 22 11  and  v(T,S)=v(Q,S) =tan  1 
(1/2 )  /rr=0 .14  76. 

If  S  and  T  are  fuzzy  sets,  call  them  ^  and  v,  we  define 
re(P,u,v)=l  if  v  (R)  2ii(P)  at  some  point  R  on  the  ray  emanating 
from  P  in  direction  0.  (Recall  [21  that  v  surrounds  ^  if, 
for  any  point  P  and  any  path  n  from  P  to  B,  there  exists  Re-; 


such  that  v  (R)>jj  (P)  . )  We  then  define  v(P,u,v)  as  ^ 

/2TTrfi  (P  ,  u  ,  v)  d0  ,  and  we  define  v(y,v)  by  taking  the  min  over 
0  0 

all  P  in  the  plane.  (In  the  case  of  taking  the  "average" 
value  of  v(P,u,v)  the  denominator  for  the  average  is 
Evidently  this  generalizes  the  crisp  definition. 

Proposition  3.1.  If  T^T' ,  then  v  (P,T)>v (P,T' )  for  any  P 
and  v  (S,  T)  >v  (S,T '  )  for  any  S. 

Proof :  This  follows  immediately  from  the  fact  that  r(P,T)~ 

r0(P,T')  for  any  P.  Analogously,  in  the  fuzzy  case,  if  v2v  '  , 
then  v  ( P ,  u  ,  v )  >.v  ( P ,  u ,  v  ' )  for  any  P  and  u  .  || 

It  is  not  hard  to  see  that  v(P,T)  is  a  continuous  function  of 
the  position  of  P.  On  the  other  hand,  v(P,T)  need  not  increase 
as  P  moves  closer  to  T,  even  if  T  is  convex,  as  illustrated  in 
Figure  10. 

If  T  subtends  angle  a  from  P,  we  evidently  have  v(P,T)= 
ci/2Tr;  thus  as  P  approaches  T,  if  T  is  convex  v(P,T)  approaches 
1/2  since  a  approaches  tt  .  Note  that  when  P^T  we  have  v(P,T)=l. 


3.2.  Topological  surroundedness 

Even  if  v (P , T) =1 ,  T  may  not  surround  P  in  the  usual  sense, 
since  there  may  be  a  curved  path  from  P  to  B  (the  "background" 
region,  outside  the  picture)  that  does  not  intersect  T,  as  illus¬ 
trated  in  Figure  11.  In  this  section  we  introduce  an  alternative 
definition  of  quantitative  surroundedness  that  is  more  closely 
related  to  the  usual  (topological)  definition. 

Intuitively,  the  degree  to  which  T  (topologically)  surrounds 
P  is  related  to  how  much  a  path  from  P  must  change  direction  in 
order  to  reach  B  without  intersecting  T.  For  example,  if  T  is 
a  spiral,  and  P  is  "surrounded"  by  a  very  large  number  of  turns 
of  T,  a  path  from  P  that  does  not  intersect  T  must  turn  through 
a  very  large  multiple  of  2r  before  it  can  reach  B. 

Let  be  any  rectifiable  path  from  P'  to  B,  starting 
at  P'  in  direction  0,  that  does  not  intersect  T.  (If  no  such 

path  exists,  we  define  t  (P  '  ,  T)  =<*>. )  Let  C  (P',T  )=/ 

0  8 

|c  (P) I dP ,  where  c  (P)  denotes  the  absolute  curvature  of 

e  11 0 

Tig  at  an  arbitrary  point  P  on  tt0.  Let  C0(P',T)=inf  (P',T); 

11  e  Q 

if  T  is  a  "well  behaved"  set  (e.g.,  a  region)  and  P ' /T  and  P’ 

is  not  inside  a  hole  of  T,  Cq(P’,T)  will  be  finite.  Finally, 

.  2tt 

let  t(P',T)=  !q  C0(P',T)d0,  the  average  of  C0(P',T)  over  v. 
[Alternatively,  we  could  have  used  110^  C0(P',T),  but  using 
the  average  allows  our  definition  to  be  sensitive  to  "partial” 
surroundedness  of  P'  by  T.  For  example,  if  T  is  a  circle  with  a 
small  gap,  and  P'  is  at  its  center,  there  exists  a  direction  in 
which  does  not  have  to  turn  at  all,  so  that  the  inf  definition 


gives  0,  as  if  T  were  not  there  at  all.  On  the  other  hand,  the 
averaging  definition  reflects  the  fact  that  some  oaths  may  have 


to  turn  by  as  much  as  it  before  they  can  get  out  of  T;  in  fact, 
the  average  is  approximately  f,  but  it  gets  smaller  as  the  gap 
in  the  circle  gats  wider.] 

If  S  is  a  ("well  behaved")  set,  we  define  t(S,T)  as 
the  inf  of  t(P,T)  for  all  PtaS.  If  S  and  T  are  fuzzy  (call 
them  n-  and  v),  we  use  exactly  analogous  definitions,  except 
that  ttq  is  a  path  from  P  to  B  such  that  v(R)<p(P)  for  all  R  on 

n  ;  this  is  the  fuzzy  version  of  "does  not  intersect  T".  In 

y 

the  fuzzy  case,  t(y,v)  would  be  ^  ^  j  ^  ^  if  we  use  the 

averaging  definition. 

Proposition  3.2.  If  T£T' ,  then  t (P,T) 2t (P,T 1 )  for  any  P, 
and  t  (S ,T ) >t  (S,T ' )  for  any  S. 

Proof :  This  follows  immediately  from  the  fact  that  any  path 

(from  any  P)  that  meets  T*  also  meets  T.  Analogously,  in 
the  fuzzy  case,  if  v2v',  then  t  (u  ,  v)  2t  (u  ,  v  '  )  for  any  u.  '! 

It  is  not  hard  to  see  that  t(P,T)  is  a  continuous  function  of 
the  position  of  P.  However,  t(P,T)  need  not  increase  as  P  moves 
closer  to  T,  even  if  T  is  convex,  as  we  see  from  the  examples  in 
Figure  10. 

Let  T  be  convex  and  subtend  angle  a  at  P.  It  is  easily 
seen  that  for  all  i  outside  that  angular  sector,  paths  from  P 
to  B  exist  that  do  not  turn  at  all  and  do  not  meet  T;  but  if 
is  inside  the  sector,  say  3  away  from  the  nearer  boundary  of 
the  sector,  a  path  from  P  in  direction  0  must  turn  by  at  least 
in  order  to  reach  B  without  meeting  T.  Moreover,  such  paths 
exist  that  do  not  turn  by  more  than  i. 


It  follows  that  t(P,T) 


is  just  the  average  value  of  S  for  all  directions  9  in  the 

sector;  this  is  evidently  just  a/2.  In  particular,  as  P  ap¬ 
proaches  T,  t  (P , T )  approaches  tt/2,  since  a  approaches  n.  Note  * 

that  when  P€T  we  have  t(P,T)=*>. 


3.3.  Surroundedness  for  digital  polygons 
|  For  subsets  of  digital  pictures  our  approaches  to  quan¬ 

titative  surroundedness  have  to  be  "digitized." 

In  the  case  of  visual  surroundedness  we  define  v ,  .  (S,T)  = 

dig 

I  min  v..  (P,T)  for  digital  polygons  S  and  T,  where  v  ,  .  (P,T) 

PtBP.(S)  019  dig 

v(P,T)=ci/2tt  if  T  subtends  angle  a  from  P.  The  rectangles  S 

and  T  in  Figure  9  may  be  considered  to  be  digital  polygons, 
l4  for  example.  Then,  the  values  of  (S,T)  and  v^.^(T,S)  re¬ 

main  the  same  as  given  in  Section  3.1  by  v(S,T)  and  v(T,S), 
respectively.  The  straightforward  approach  to  computing 

*  v  ,  .  (S , T)  would  be : 

dig 

angle  =  +<», 

compute  the  convex  hulls  S'  ,T'  of  S,T  using  any 
desired  linear  time  algorithm 
for  all  points  P  in  BP(S’)  do 

1.  compute  the  two  tangents  from  P  to  T ’  and 
the  angle  a  between  these  tangents, 

^  2 .  Lf  angle  >a  then  angle  =  a  od 

return  angle/2^. 

Since  when  P  moves  around  S'  the  related  tangential  points 

>  Ql ,Q^ fBP (T 1 )  move  around  T*  monotonical ly ,  this  straightforward 
algorithm  leads  to  an  0(N)  time  algorithm,  >I=card  BP(S)  + 

card  BP (T ) ,  by  using  two  points  to  the  actual  tangential  points 

>  in  BP(T'). 

In  the  case  of  topological  surroundedness,  besides  the 


restriction  of  aS  to  BP(S)  we  have  to  digitize  the  set  of 
possible  directions  -  for  paths  from  BP(S)  to  the  background  B, 
for  a  digital  polygon  S.  Let  us  assume  that  is  restricted 


to  the  set  ang  ={n*  — :  n=0 , 1 , 2 , . . . ,m-l } ,  for  mil.  Then 
hn  m 

C  ,  { P  ,  T )  denotes  the  minimal  angle  that  a  path  n  in  direc¬ 
ts  6 

tion  e  starting  at  P  may  take  around  T  to  B ,  as  defined 

in  Section  3.2,  and  t™ .  (P,T)  is  defined  by  —  Z  Cr  (P,T)  , 

dig  m  y 

^  6  t ang 

m 

for  a  digital  polygon  T  and  a  point  P.  Finally,  we  have 

t1? .  (S,T)  =  min  t™ .  (P,T).  Obviously,  the  computational 

dlg  P£BP (S )  dlg 

requirements  for  computing  t™^  (S,T)  exceed  those  for  com¬ 
puting  the  visual  surr oundednecs  measure  (S,T),  but 

still  t™^  seems  to  be  a  practically  useful  function.  For 

O 

example,  in  the  situation  of  Figure  9  we  have  t^  (P,T)  = 
-|(tan-1  (5/6)+0+0+0+0+0+0+0)=0.0276-n  =  0.0868.  It  follows 
that  t^ig  (S,T) -t^ig(P/T)=0. 0868.  Analogously,  t^ig(T,S)s 

Q  „ 

t^g  (Q,  S)  =0 . 0184  •  n  =  0 . 05  79  .  For  computing  t^  (S,T) 
nearly  the  same  algorithm  may  be  used  as  for  computing 
vd^g(S,T),  but  with  some  extensions.  After  computing  the  two 
tangents  from  P  to  T '  we  determine  a=t™^g(P,T')  by  using  the 
minimal  angular  differences  to  these  tangents  if  -3  is  between 


these  tangents;  otherwise  C.(P,T')=0.  Thus,  t^  (S,T)  with 
N=card  BP(S)+card  BP(T)  may  be  computed  within  0(mN)  time  in 


the  worst  case  sense. 


We  have  given  definitions  for  quantitative  adjacency  and 
surroundedness  for  crisp  and  fuzzy  sets  in  the  real  plane,  and 
also  for  crisp  digital  polygons;  we  have  not  considered  the 
digital  fuzzy  case.  We  have  briefly  described  some  algori¬ 
thms  for  computing  “these  quantities  in  the  digital  case,  but 
many  issues  have  been  left  open;  for  example,  fast  algorithms 
for  the  adjacency  measure  in  the  general  case  (arbitrary 
polygons)  are  still  an  open  problem. 

The  proposed  measures  should  be  of  interest  in  the 
study  of  stochastic  geometry  in  the  real  plane.  At  the  same 
time,  these  measures  can  be  used  to  characterize  relation¬ 
ships  between  objects  in  a  segmented  digital  picture,  or  to 
compare  objects  in  two  different  pictures. 
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Figure  4.  The  degree  of  adjacency  of  a  region  T  to  a 
point  P  is  not  necessarily  a  monotonically 
decreasing  function  of  d(P,T),  and  is  not 
necessarily  a  continuous  function  of  the 
position  of  P. 
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Figure  5.  Counterexample  to  the  fuzzy  generalization 
of  Proposition  3. 


END 

FILMED 

8-84 

DTIC 


